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Coherent Synchrotron Radiation can severely limit the performance of accelerators designed for 
high brightness and short bunch length. Examples include light sources based on ERLs or FELs, and 
bunch compressors for linear colliders. In order to better simulate Coherent Synchrotron Radiation, 
the established 1-dimensional formalism is extended to work at lower energies, at shorter bunch 
lengths, and for an arbitrary configuration of multiple bends. Wide vacuum chambers are simulated 
by means of vertical image charges. This formalism has been implemented in the general beam 
dynamics code Bmad and its results are here compared to analytical approximations, to numerical 
solutions of the Maxwell equations, and to the simulation code elegant. 



I. INTRODUCTION 

It is envisioned that future accelerators will call for 
shorter beams of higher intensity. A possible limiting 
factor in these efforts is an increase in energy spread and 
transverse emittance, as well as a micro-bunching insta- 
bility, due to Coherent Synchrotron Radiation (CSR). 

The first CSR calculations were performed by 
Schwinger in 1945. Using a Green's function method, he 
arrives at the power spectrum of a single charge bend- 
ing in free space as well as between infinite conducting 
plates, and thereby computes the coherent power radi- 
ated by a collection of charges [1]. Warnock extends this 
work to include the longitudinal impedance on a bunched 
beam 0. Many papers covering the history and impor- 
tance of CSR forces can be found in 0|. 

This paper uses an approach to calculate the CSR 
wake-field originating with Saldin et al. [4| and gener- 
alized by Sagan Here we calculate the CSR force 
between two charges traveling on the same curve, and 
integrate over a longitudinal bunch distribution to give a 
longitudinal wake-field. Transverse particle coordinates 
and transverse force components are neglected. The for- 
malism developed here is generalized to include arbitrary 
lattice configurations of bends and drifts, including, for 
example, radiation from one bend entering another, and 
bend radiation extending into drift regions. 

Simulating CSR effects is the subject of a number of 
codes. The method here is implemented in the particle 
tracking code Bmad [6]. Our simulation results are com- 
pared with two of the codes described in [7[ and with 
approximate analytic formulas. 
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FIG. 1: A particle at point P' kicks a particle at point P. 



to the source particle at point P' and retarded time t' 



E(P) = 



4(L - W) + ^L x [(L - W) x a'] 
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It will be assumed that both particles have the same con- 
stant speed /3 = v/c, n' and n are the unit velocity vec- 
tors for the source and kicked particles respectively, and 
L = P — P' is the vector from the source point to the kick 
point. The retarded time t' is related to t via t — t' — L/c. 
At time t, the source particle has a longitudinal position 
z' with respect to the bunch center and the longitudinal 
position of the kicked particle is z. The distance £ = z— z' 
between the particles at constant time can be computed 
via the equation 



C = L S -(3L 



(2) 



II. TWO PARTICLE INTERACTION 

In order to compare the force acting on one particle 
from the radiation that is emitted by another, the analy- 
sis starts by considering two particles of charge e follow- 
ing the same trajectory as shown in Fig.[TJ The Lienard- 
Wiechert formula [8] gives the electric field E(P) at the 
position of the kicked particle at point P and time t due 



where L s is the path length from P' to P. Generally, the 
relativist ic approximation (3 = 1 will be made. However, 
some terms in 1 — (3 ~ 1/2 7 2 will need to be retained. 

The first term on the right hand side of Eq. (pQ) has 
a l/( 2 singularity at small distances. Following Saldin 
et al. [J], this singularity is dealt with by dividing the 
electric field into two parts. The space charge compo- 
nent E sc , which contains the singularity, is the field that 
would result if the particles where moving without ac- 
celeration along a straight line. The CSR term, E CSR , is 
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FIG. 2: i^csR (left) and (right) as a function of £ for a bend. 



what is left after subtracting off the space charge term 
e n 
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Ecsr = E — E s 



(3) 



The rate K = dS/ds at which the kicked particle is 
changing energy due to the field of the source particle 
is 

K = if CSR + if sc = e n • E CSR + e n • E sc . (4) 

Following Saldin et al. [4], the transverse extent of the 
beam will be ignored in the calculation of K CSR . How- 
ever, the inclusion of the finite beam size will be needed 
to remove the singularity in the calculation of K sc as 
discussed in Section HVl 



III. CSR CALCULATION 

The source point P' and the kick point P will, in gen- 
eral, not be within the same lattice element. Because 
the transverse extent of the beam is being ignored, all 
elements will be considered to be either bends or drifts. 

In Fig. [TJ R is the bending radius and g = 1/R is the 
bending strength of the element that contains the source 
point P'. The magnitude of the acceleration is a' ~ c 2 /R. 
This element ends at point O. The angle between P' and 
O is 0, and d = R<p is the path length between P' and 
O. 

Between point O and the kick point P, di is the path 
length within the i th element, i = 1,...,7V, where TV 
is the number of elements in this region. For the last 
element, d^ is the distance from the start of the element 
to point P. For the i th element, fa is the bend angle, Ri 
is the bend radius, and gi = 1/Ri is the bend strength. 
For a drift fa = gi = 0. 

In Fig. [TJ (v,w) are the coordinates of point P with 
respect to point O with the v-axis parallel to the orbit's 
longitudinal s-axis at point O and the w-axis pointing 
upwards towards the inside of the element containing the 
point P'. 



With this notation, the difference in v and w from the 
beginning of an element to the end is 



Awi 



Ri (sm(fa + ipi) — sin ipi) for a bend 

di cos ipi for a drift 

Ri (cos ipi — cos(fa + ijii)) for a bend 

di sin ipi for a drift 



(5) 



where ipi is the orientation angle at the entrance end of 
the element 



i-l 

k=i 



(6) 



The above formulas are able to handle negative bends 
(beam rotating clockwise). For a negative bend Ri, gi 
and fa are negative while di = Ri fa is always positive. 

With the assumption that all bend angles are small, v 
and w can be approximated by 



where 



v = v\ - V3 , 
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W = UJ>2 , 
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and the small angles have been retained to second order. 

The angle of the vector n with respect to the v-axis 
is 9 = YliLi 9i di. In terms of v and w : the components 
of the vector L are 



L v = v + R sin (j) = [v\ + d] 



^3 + ~ 

gd 2 

L w = w — R (1 — cos fa) = U02 — , 



i 2 d 3 



(9) 



L = y/L*+Ll 
= Wi + d] - 
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Again angles are retained to second order. The path 
length is simply 



N 



= d + ^2 di = d + ^i . 



This, With Eq. ©, gives 



c 



2 7 2 



^3 



^ 2 d 3 1 (2u 2 -gd 2 ) 2 
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(10) 



(11) 



where terms to second order in combinations of angles 
and I/7 are retained. Substituting these expressions into 
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Eq. (pQ), and defining 
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1 



u 2 +gdvi + -gd' 



(12) 



k = 7 (0 + d) , r = 7 (d + z/ x ) 

the individual terms in Eq. (pQ) read as 

1 87 9 r 3 
(L-L-/?n') 3 " (r 2 +a 2 ) 3 



(13) 



n- (L-L/3n / )/7 2 = 
n • (L x [(L - L/3ri) x a 7 ]) /c 2 = 
Putting these together yields 

(r 2 - a 2 ) (a — TK) 
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2 7 5 r 

# (r 2 — a 2 ) (a - r k) 
27 5 r 
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(r 2 +a 2 ) 3 
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(r 2 +a 2 ) 3 



(14) 
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While we have used S.I. units, the classical radius r c and 
the mass m has been used to make the formula indepen- 
dent of the unit system. 

From Eq. ([HP , K° SK , which is K CSK restricted to the 
special case where points P and P' are within the same 
bend, reduces to Eq. (32) of Saldin et al. [4|], 



4 r c mc 2 7 4 
~R 2 



^ 2 /4-l 
[2(l + 2 /4) 3 

•3 2 /4 1 



(15) 



(l + </> 2 /4) 3 (1 + 2 /12) 2 



where = 70, a = R(j) 2 /2, n = </>, and r = Rcj). This 
equation is valid for > 0; for < 0, if C sR is, to a very 
good approximation, zero. 

In the limit of small £, ^c S r nas a limiting value of 



2 ^4 



^csh(C) 



-4r c rac 7' 
3i? 2 



for C < -o 



(16) 



At large values of £, ^ is cubic in </> so that 
(24C/i^) 1/3 . With this, Eq. becomes 



^csr(C) 



2 r c mc 2 

3 4/3 #2/3 £4/3 



for C > -o , 



(17) 



which corresponds to Eq. (10) of Saldin et al. [4] (note the 
error in the denominator of Saldin et al. Eq. (10)). Fig- 
ure [2] shows if ° SR (£), which changes sign at ( « 1.8 R/j 3 . 
The long tail at £ > 1.8i?/7 3 cannot be neglected since 
the integral J °° dC-Kcsii(C) is zero. The vanishing of the 
integral is a reflection of the fact that a closed loop of 
charged particles of uniform density does not radiate. 
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FIG. 3: /csr as a function of f for a bend. The dashed line 
is the large f approximation as given in Eq. (|21[) . 



The fact that K CSK is highly peaked in amplitude 
near £ = can be problematic for simulations at ultra- 
relativistic energies because the characteristic longitudi- 
nal distance between particles or mesh points needs to 
be less than R/j 3 . One way of dealing with the peaked 
nature of K CSK is to first consider the kick from a line of 
particles of density X(z) and then to integrate by parts 



f-) = rdz'\(z')K CSR (z-z') 

J CSR ^-OO 

■/ 



dz ,d>^)_ _ z , 
dz' 



(18) 
(19) 



where 



W* -*') = - / Z dz " K OSR (z - z") 

J — CO 



(20) 



7° SR , which is 7 CSR for P and P' in the same bend, 
is plotted in Fig. 03 The peaked nature of K CSK has 
been smoothed over at the cost of having to deal with 
a derivative of A. For £ ^> R/j 3 , the approximation 
of Eq. (fTTj) can be used to calculate an explicit ultra- 

1, 



relativistic equation for 7° SR as in 



ro (n = -2r c mc 2 1 

1 csn\^) 31/3^2/3 £1/3 



r » R 

for C > -o 



(21) 



Equation ([2T]) is also plotted in Fig. [3l 

While, in general, it is helpful to have explicit formu- 
las, for the purposes of evaluation within a simulation 
program this is not needed. The alternative is to use an 
exact implicit solution. Because Eq. ([TTj) and Eq. (J4j) 
are rational functions, Eq. ([20]) can be integrated. The 
last term that compensates K sc can be integrated to 
r c rac 2 /(7 2 C). The other terms can be written with 
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FIG. 4: Comparison between Eq. (|28|) and an exact integra- 
tion of Eq. (I2S)) . 
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With dr/dd = 7, da/dd = jgr, and dn/dd = 75, one 
can further simplify to 



2r c mc 7 I — J — 



dd 
d_ 



d 
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This can be integrated over £ to yield 
'2 7 (r + a ft) 



JcsrO,^) -r c mc 2 



a* 



1 

7 T C 



Tor • (25) 



It can be shown that, while quantities like d and g are 
discontinuous across element boundaries, r, a, and ft are 
continuous and hence 7 CSR is a continuous function as it 
should be. 

Equation ([25j) is the main result of this paper. Us- 
ing Eq. ([25]) , the integration of Eq. (fT9|) in a simula- 
tion program can be done via interpolation of Eq. (fTT]h 
Equation ([25]) has several advantages over equations like 
Eq. ([2T]) . It is applicable at lower values of 7 s that is, 
at lower energies and/or smaller length scales. Addition- 
ally, Eq. (|25|) has no singularity at small £, and it can be 
used to handle any combination of elements between the 
source and kick points. 



IV. SPACE CHARGE CALCULATION 

The singularity at small £ in the space charge term 
E sc is removed by considering the finite transverse beam 
size. This term is equivalent to the problem of calculating 



the field given a static distribution of charges. It will be 
assumed that at any longitudinal position the transverse 
profile of the beam is Gaussian. Thus, a longitudinal slice 
of the beam will produce an energy change for a particle 
at longitudinal z and transverse offset (x, y) from the slice 
center of 

/CO pOO 
I dx' dy' p(x' ,y' ,z')dz' (26) 
-00 J — 00 

r c mc 2 7 z 



( 7 2 (z - z') 2 + (x - x') 2 + (y - y'yy/z 



where p is the bi-Gaussian distribution 



p(x,y,z) 



P(*) 



2na x a y 



exp 



x 



2o* 



(27) 



A heuristic solution for Eq. 

(x < 3 a x and y < 3 a y ) is 



in the region of interest 



dK s 



r c mc 2 sign(C) p{z')dz' 



a x (j y exp 



7lCI+7 2 C 2 



y 2 , 

(28) 

-1 for negative z, and 



where sign(z) is 1 for positive and 
C = z-z'. 

Equation ([28]) is exact in the limit z = and z — > 00, 
and is an excellent approximation in the region in be- 
tween. This is illustrated in Fig. [4j which shows K sc 
as a function of z as computed from an integration of 
Eq. ([26|) and from the approximate Eq. ([28)) . The par- 
ticular parameters chosen for the computation are given 
in the figure. Two cases were considered. One where 
the kicked particle is on- axis, and the other where the 
kicked particle is displaced by cr x off-axis. As can be 
seen, Eq. ([28j) gives an excellent approximation to the 
longitudinal space charge kick. 

At high energies, the CSR energy kick is independent of 
the beam energy as indicated by Eq. ([2T|) . On the other 
hand, the factors of 7 in the denominator of Eq. ([28]) 
insures the K sc will be a decreasing function of 7. Past 
some point in an accelerator, the increasing 7 will make 
the effect of K sc small compared to the effect of if C sR- 
To estimate that point, consider the maximum energy 
kick in a Gaussian bunch of N e particles with 



p(z) = 



'2 tt a. 



exp 



z 



(29) 



The space charge kick is maximum at x, y = and using 
Eq. (|2SD gives 



i^sc(0,0,z) 



N e r c mc 2 



(30) 
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7C + 7 2 C 2 



The dominant term in the denominator in the inte- 
grand is either a x a y for small or j 2 ( 2 for large (. 
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As an approximation, the middle term in the denom- 
inator will therefore be ignored. The approximation 
sinh (z^/a 2 ) ~ z C/ cr z wml a l so be made. This approx- 
imation is justified since, as shown below, the region of 
interest for z is around z « a z and for realistic beam 
parameters, a x a y <C 7 2 o" 2 an d therefore the only signifi- 
cant contribution to the integral will come in the region 
C < cf z . With these approximations, Eq. ([3Q|) becomes 



K sc (z) 



2 N e r c mc 2 z exp 



2<r 



f 

Jo 



C exp 



2(7 



(31) 

The integral can be evaluated to (2j 2 ) 1 e a T(a) where T 
is the exponential integral and a = a x cr z /(2j 2 al). For 
small a this is approximately — ln(a). The maximum 
if sc ,max occurs at z = a z . Using a x a y <C 7 2 cr 2 gives 



if sc ,max 



The maximum CSR kick at large energies is, from 
Eq. CHj), 




ifcsR, max « J 



iV e r c m c 

(iZM) 1/a 



(33) 



The ratio is 

if gc ,max 
i^ CSR ,max 



0.3 log 



r 



2/3 



(34) 



The condition for K sc begin small can be written as 

1/2 / d \ 1/3 



7 > 



0.3 log 



2 7 V 



(35) 



The condition for this can be well approximated by not- 
ing that the quantity in square brackets on the LHS of 
Eq. ([35]) is slowly varying and never extremely large. 
Thus the required condition is 



nr I R 

7 > M — 
\cr z 



1/3 



(36) 



where M is a number of order unity. In this case, <j z must 
be interpreted as the characteristic longitudinal distance 
over which the bunch density is changing. 



V. CSR IN BMAD 

The above algorithm for simulating CSR and Longitu- 
dinal Space Charge (SC) has been implemented as part 
of the Bmad [6] subroutine library for relativistic charged- 
particle simulations. Bmad simulates a beam as a set of 
particles. The beam is tracked through a lattice element 
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FIG. 5: Bmad implementation of the CSR algorithm. The 
beam of particles is divided up into a number of bins. The 
contribution of a particle to a bin's total charge is determined 
by the overlap of the particle's triangular charge distribution 
and the bin. 



by dividing the element into a number of slices. Tracking 
through a slice involves first propagating the particles in- 
dependently from each other and then applying the CSR 
and SC energy kicks. To calculate the energy kick, the 
beam is divided longitudinally into N b bins as shown in 
Fig. [5l For computing the charge in each bin, each beam 
particle is considered to have a triangular charge distri- 
bution. The overlap of the triangular charge distribution 
with a bin determines that particle's contribution to the 
total charge in that bin. The width of the particle's tri- 
angular charge distribution and the number of bins are 
set by the user. The bin width is dynamically adjusted 
at each time step so that the bins will span the bunch 
length. Increasing the particle width smooths the distri- 
bution at the cost of resolution. 

The charge density Xi at the center of the i th bin is 
taken to be A$ = pi/Az b where pi is the total charge 
within the bin and Az b is the bin width. The charge 
density is assumed to vary linearly in between the bin 
centers. The CSR energy kick for a particle at the center 
of the j th bin after traveling a distance ds slice according 
to Eq. (HSJ is then 



d£j = ds sUce (Xi — Xi-i 

i=l 

where 

-^CSrO) = I CI 



IcsrU - i) + IcskU - i + 1) 



t (z = j Az b ) 



(37) 



(38) 



Evaluation of I C sn,j involves inversion of Eq. ([TTj) to 
obtain d. Because z is a monotonic function of d, New- 
ton's method [To[ is used to find numbers d\ and di which 
bracket the root and then Ridders' Method [10[ is used 
to quickly find d. 

In deriving Eq. (j37j), the approximation 

^ + ^ b J T ( \ ^ A -^CSr(j') + ^CSr(J + 1) / on x 

dzI CSK (z) w Az b (39) 

j Az b Z 

has been used. Generally this is an excellent approxima- 
tion, except when j = and Az b ^> R/j 3 , as shown in 
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Fig. 03 Here, however, the integral can be done exactly 
assuming that the source and kick points lie within the 
same element 

Jo dzI - {z) = 7^ n \J(K^) —■ 

(40) 

Once the energy kick at the centers of the bins is calcu- 
lated, the energy kick applied to a particle is calculated 
via interpolation assuming a linear variation of the kick 
between bin centers. 

In calculating the energy kick, The computational time 
for calculating the charge in the bins charge scales as N p , 
the number of particles in the simulation. The compu- 
tational time for calculating the energy kick at the bin 
centers scales as N 2 , and the time for calculating the 
energy kick of the particles scales as N p . 



where E^ is a complex 2-dimensional vector related to 
the perpendicular electric field, po is the charge density, k 
is the wave number, and R is the magnet bending radius. 
It is solved using a finite-differencing method. 



VII. CSR IN ELEGANT 

The particle tracking code elegant (version 17.2.2) 
uses Eq. (|^3|) to compute CSR kicks without shielding 
by a vacuum chamber [l2|. The charge distribution X(z) 
and its derivative dX/dz are calculated by binning the 
macroparticles and then employing a smoothing filter. 



VIII. CSR WAKE FORMULA 



Chamber Walls 



A. Transient Effects at Magnet Edges 



The simulation incorporates the shielding of the top 
and bottom chamber walls by using image currents. Ap- 
pendix El explains why neglecting the width of a chamber 
is a good approximation when it's width is larger than its 
hight. Here, because the image current is well separated 
from the actual beam, there are no singularities to deal 
with, K sc does not have to be subtracted, and a straight 
forward integration is done using Eq. ([T8]h and Eq. (jHj). 

^•( imag e) =2*ds slice (41) 

Ni N b 

x £ (-!)* £ ft * K ( z = 0" - i) Az b, V = kh) , 

k=l 2=1 



where qi — \ Az^ is the charge in a bin, h is the cham- 
ber height, and k indexes the image currents at vertical 
displacement y = ±kh. The number of image charges 
Ni needs to be chosen large enough so that the neglected 
image currents do not have a significant effect on the sim- 
ulation results. Because the relevant angles are not small, 
the image charge kick K must be calculated without the 
small angle approximation, as in Eq. (Q3J). 



VI. AGOH AND YOKOYA CSR CALCULATION 

Agoh and Yokoya (A&Y) have developed a code to cal- 
culate CSR wake fields by directly integrating Maxwell's 
equations on a mesh representing a rectangular beam 
chamber [11]. The approach depends on the paraxial 
approximation, a rigid Gaussian bunch density, small 
chamber dimensions relative to the bending radius, and 
ultra-relativistic particles. Using a co- moving coordinate 
system (x,y,s) in Fourier space, they are able to reduce 
the problem to a tractable two-dimensional differential 
equation, 



d „ i 
— Ei = — 

ds 2k 





2k 2 x\ 


"(vi- 


f * ) 





E ± V ± p 



, (42) 



Using retarded fields, Saldin et al [4[ derive, in the 
ultra-relativistic limit, a formula for the wake-field due 
to a bunch entering from a drift region into a bend 



dS 

ds 



2N e r c mc 2 J X(s - s L ) - X(s - 4s L ) 



1/3 



Jss L {s-s<y/*ds> as } ' 



(43) 



SL 



Rf 
24 



where (j) is the angle traveled into the magnet by the 
bunch center. Eq. (|^3|) reduces to the free space steady- 
state wake-field of Eq. ([2T]) in the limit sl — » oo. 

As worked out by Emma and Stupakov @ , synchrotron 
radiation will continue to propagate and affect the bunch 
beyond the end of a bending magnet. For a finite magnet 
of length L m , an ultra-relativistic bunch at a distance x 
from the end of this magnet experiences the free space 
exit wake-field 



^ | = N e r c mc z { 4 



2 f a X ( s ~ As(L m )) 



(44) 



Jo 



As(l) 



24R 2 l + x ' 



(45) 



B. Steady State CSR in a Bend 

CSR effects in a vacuum chamber have been computed 
by the Green's function of grounded parallel plates p], 
[2[. These formulas are difficult to compute numerically, 
due to the presence of high order Bessel functions, so we 
will use an excellent approximation developed by A&Y 
PH . The impedance for the steady-state in a dipole with 
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TABLE I: Parameters used in simulations 
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K{m) 


L(m) 
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A 
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1.0 


34.0 


28.0 
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2.48 


C 


10.0 


1.0 


10.0 


10.0 


0.3 


2.48 


D 


10.0 


1.0 


60.0 


10.0 


0.3 


2.48 


E 


1.20 


0.419 


0.762 


0.3 


0.036 


3704.1 


F 


2.22 


0.678 


1.71 


2.54 


1.0 


29.22 


G 


87.9 


6.574 


8.0 


4.0 


0.3 


12.52 



horizontal plates separated by a distance h is 

1/3 00 
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+ xAi(x)(Ai(x) -iBi(x)) , 



(46) 
(47) 

(48) 



where Zo = c/io is the free space impedance, fc is the wave 
number, and Ai and Bi are Airy functions. The parallel 
plate wake-field due to a bunch with longitudinal density 
X(z) is obtained by Fourier transform: 



d£ 

ds 



pp 



-N e r c mc 2 )i( - y Z(k)X(k)e iks dk 



/CO 
X(z)e- ikz dz . 
-CO 



(49) 
(50) 



IX. COMPARISON BETWEEN BMAD, THE 
AGOH AND YOKOYA, AND ELEGANT 

In order to validate our method, we compare simu- 
lations from Bmad to those using the A&Y code and 
elegant. For ease of reading, all magnet and bunch 
parameters used are enumerated in Tab. U and will be 
referred to by a letter. 

Note that all electric fields in the graphs are normalized 

by 



En = 



2N e r c mc 2 



27r(3i? 2 <j4)i/3 



(51) 



which approximately describes the maximum amplitude 
of the CSR-Wake d£/ds. 

All simulations use a bunch charge of InC, and an 
energy of 1 GeV, unless otherwise noted. We used Bmad 
with the following parameters: number of bins Nb = 800, 
number of macro-particles N p = 4 • 10 5 , number of image 



charges Ni = 32, and tracking step ds s 



1mm. The 



triangular bin width, as in Fig. is 32 bins. Figure [6] 
shows that this choice of and N p is reasonable by 
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FIG. 6: The number of particles and the number of bins are 
varied in Bmad using parameter set A. 



varying N p and A/5, using parameter set A. While reduc- 
ing the number of particles leads to a visibly less smooth 
wake field. The number of bins has been increased, un- 
til the wakefield starts becoming less smooth. The other 
parameters were similarly varied to find the applied pa- 
rameters. 



A. Steady State Case 

Figure [3 shows the steady-state CSR kick in a bend as 
a function of z for various values of the chamber height. 
The parameters used correspond to set A of Table U and 
are the same as used for Fig. 1 of A&Y [111 ]. 

Figure [71 shows excellent agreement among the 
CSRmesh code, Bmad, and the CSR-Wake formula 
Eq. (|49|) . Note that Bmad computes the wake by tracking 
a bunch, and therefore shows the result only in the range 
of the length of the bunch. This is not problematic, be- 
cause the wake only influences particles in the region of 
the bunch. 



B. Transient Case 

Parameter sets B, C, and D are used to explore the 
transient case where some of the kick is generated from 
particles in the drift region before the bend. Set B cor- 
responds to values used in Fig. 3 of Agoh and Yokoya 
pj| , which has a chamber of relatively large cross section. 
Figure [8] (left) shows agreement between Bmad, the A&Y 
code, and the CSR-Wake formula Eq. (|43]) . However, for 
parameter set C in FigJS] (right), where the chamber is 
smaller as in Fig. 4 of [11], a discrepancy appears in the 
front of the bunch. This discrepancy is noted in [ll| and 
explained to be due to the finite chamber size. 

We use parameter set D to analyze whether the cham- 
ber width is responsible for this difference, and one sees in 
Fig. [9] that the effect remains and thus appears to be due 
to the reduced chamber height. It also is larger for longer 
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FIG. 8: The transient case for a large chamber using parameter set B with a large chamber (left), and set C with a smaller 
chamber (right). The A&Y code (dots), Bmad (circles), and the CSR-Wake formula Eq. ([43]) (lines) agree well in the former 
case, but differ in the latter case. 



magnets. Varying the number of Fourier coefficients used 
in the A&Y calculation in Fig. [10] does not change this 
result either, verifying that we chose a reasonable number 
of Fourier coefficients for solving the Maxwell equations 
numerically. Reversing the set D height and width, as in 
Fig. HU somewhat reduces the discrepancy. Thus, again 
indicating that the effect is apparently due to the reduced 
chamber height. Practically speaking, the difference in 
these wake-fields only becomes appreciable far in front 
of the bunch, where there are few particles to affect. To 
see how strong the deviation becomes for especially small 
chamber heights, we repeat parameter sets B and C with 
h = 2cm in Fig. [12] and the disagreement between Bmad 
and the A&Y code is again not very large but significant. 



C. Realistic Magnets 

To evaluate how significant the differences are in real- 
istic magnets, we use parameter sets E, F, and G which 
correspond to the JLab TH2, CESR Analyzer magnet, 
and Cornell ERL CE magnets, respectively. Wake-fields 
for the first two are plotted in Fig. [T3] showing good 
agreement between the A&Y code and Bmad. The free 
space steady-state wake-field Eq. ([43]) and results from 
elegant are plotted for reference. 

The principal detrimental effects of the CSR-Wake in 
an accelerator are energy loss and increase in energy 
spread of a bunch. The transverse bunch distribution 
can also be damaged and is mostly influenced by in- 
creases of the energy spread which, through dispersive 
orbits, couples to transverse motion. To visualize CSR- 
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FIG. 9: Length dependence for the transient cases using the 
wide chamber in parameter set D. Magnet lengths of 60cm 
(top), 80cm (middle), 100cm (bottom) show that the codes 
deviate for longer magnets, but Bmad (circles) agrees with 
CSR-Wake Eq. (|43[) . Thus, reduction of the chamber height 
and longer magnet length lead to discrepancies of codes. 



driven energy loss and energy spread in a single magnet, 
we plot in Fig. [TH the average and RMS electric field 
across the bunch distribution as a function of distance 
into the magnet. This is done using parameter set G, 
where we examine how well one can ignore the chamber 
width, as done in the Bmad calculation. One sees that, 
in this shielded case, the chamber width w in the A&Y 
code begins to change the wake field when it is compa- 
rable or less than 4cm, the height of the chamber. This 
effect is heuristically explained in appendix [A[ It again 
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FIG. 10: Using parameter set D, varying the number of 
Fourier coefficients (Nk) used in the A&Y calculation does 
not change the numerical solution of the Maxwell equations, 
and therefore does not account for deviations in the presented 
comparison of codes. 



shows that ignoring the chamber width, as in Bmad is 
a reasonable approximation when the chamber is wider 
than high. 

It is apparent that shielding by a vacuum chamber 
reduces the power emitted by CSR very effectively. It 
does not reduce the energy spread nearly as much, and is 
therefore not as effective for preserving bunch properties 
as one might have concluded from the reduced radiation 
power. Interestingly, the 4cm wide chamber even pro- 
duces larger RMS E s than wider chambers. 



D. Exit Wake 

The method in this paper can correctly account for 
the CSR-Wake in a drift section following a magnet. Us- 
ing parameter set A, Fig. [15] shows this wake as a func- 
tion of the distance d from the end of the magnet in 
free space and with shielding. The free space case shows 
excellent agreement between Bmad and the CSR-Wake 
formula Eq. (EHj) . 

The other codes do not allow simulations with shield- 
ing in this regime. We therefore compare Bmad with a nu- 
merical solution using retarded fields. In order to avoid 
the complication of computing retarded times and po- 
sitions of the bunch distribution and its image charges, 
we do not use Lienard-Wiechert fields, but rather Jefi- 
menko's equations Q. In general, the electric field due 
to a 1-dimensional charge density p(s, t) and current den- 
sity J(s,t) at a position s, time t, and height h is 



E(s,t;/i) 



1 



47T60 



ds' 



(52) 



+- 



L 2 



\dp(s',f)-\ 


1 


~dJ(s',fy 










1 



where L is the vector from position s' to position s 
at height h, L is its magnitude, and the brackets [] are 
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FIG. 11: Length dependence for the transient cases using a 
tall chamber by reversing width and height in parameter set 
D. Magnet lengths are 60cm (top), 80cm (middle), and 100cm 
(bottom). 



evaluated at the retarded time t' = t — L/c. As com- 
pared with integrating over retarded fields at s' , using 
Jefimenko's equations has the advantage that one never 
has to solve for the retarded time in equations of the type 
t f = t-\r-r(t f )\/s. 

The total electric field due to alternating image charges 
is then 



E imag es(s, t) = 2 ^(-l) n E(s, t ; n h). (53) 

n=l 

Applying this to the geometry of a bend followed by 
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FIG. 12: The same as Fig. [8] except with the chamber height 
reduced to 2cm. A&Y code (dots), Bmad (circles), and the 
shielded steady-state CSR-Wake Eq. (09} (line) show slight 
deviations from one another. 



a drift, Fig. [T5l (bottom) shows excellent agreement with 
Bmad. 



E. Coherent Energy Loss 

Our final test of Bmad compares the total coherent 
energy lost for various particle energies and chamber 
heights to the integration of the power spectrum. In 
general, for N particles traveling on the same curve at 
different phases <j) ni the N particle power spectrum is 



duo 



N 



n=l 



duo 



(54) 



/ 



N 



N- 



1 + (N - 1) 



duo 



duo 

A N(N — 1) 



(55) 
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FIG. 13: Realistic magnets: Parameter set E (JLab TH2 
magnet) line (top), set F (CESR Analyzer magnet) (bottom). 
Bmad agrees with the CSR-Wake formula Eq. (|49[) better than 
the other codes at the bunch tail. 
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FIG. 15: Exit wake-field as a function of the length d from 
the end of the magnet for free space (top) and with shielding 
(bottom). Bmad (dots) shows excellent agreement with CSR- 
Wake formula Eq. (|44l (top lines) in the free space case, and 
with numerical integration over image bunches using Eq. (153 
(bottom lines). Parameters set A is used. 



resulting coherent part shows excellent agreement with 
Bmad for energies down to 5MeV and chamber heights 
down to 2mm. At smaller heights, the number of im- 
age layers used in the simulation (Ni = 64 here) is not 
sufficient to correctly model the CSR. 

In absence of shielding plates, Eq. ([55]) can be inte- 
grated exactly for a Gaussian distribution using the well- 
known free space single particle power spectrum. With a 
standard deviation <j z , the total power lost by N particles 
is 



FIG. 14: Average and RMS longitudinal electric fields for 
parameter set G (the Cornell ERL's CE magnet) for vari- 
ous chamber widths using the A&Y code, compared to Bmad, 
which has infinite chamber width. 



where dP^/dw is the single particle power spectrum, 
and X(z) is the longitudinal particle distribution. The 
first term in Eq. ([55]) is the incoherent power spectrum, 
while the second is the coherent power spectrum. In the 
presence of conducting parallel plates, the single particle 
power spectrum is given in Eq. (47) of [1]. Using this, 
Eq. ([55]) can be integrated numerically, and in Fig.[T6lthe 



P (N) = p (i) N + p (i) N ( N _i) T ( ^f.) , (56) 
where 

P W = \r c mc^ (57) 
is the power lost by a single particle, and 

rw-aT^GpWsO-is 1 - (58) 
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This formalism has been implemented in Bmad. We 
show that the longitudinal wake field compares well with 
the A&Y code of A&Y and the CSR-Wake formula 
of Warnock [2j] for the steady state, with and without 
CSR shielding by parallel plates. In the transient case 
where the A&Y code often does not follow the CSR-Wake 
formula of [13[ exactly, Bmad does agree well with that 
formula. 
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APPENDIX A: HEURISTIC SHIELDING 
ARGUMENT 



FIG. 16: Average energy losses versus shielding height for 
various energies. The bottom plot is the same as the top 
highlighting the smaller heights. Bmad (dots) agrees well with 
numerical integration of Eq. ([55]) (curved lines) using the 
shielded power spectrum of [l| down to 2mm. Steady state 
losses computed using Eq.[56]are indicated as horizontal lines. 
Parameters set A is used, with the number of image layers 
Ni 64. 



This result agrees well with Bmad in Fig. [16j The func- 
tion T(a) can be expanded asymptotically, giving the 
leading order coherent energy change in a length L as 



CSR Induced Average Energy Change (eV) 



AE (N) 



r 

-N 2 r c mc 2 — ^ 



L 



1/3- 



(59) 



This is consistent in the scaling and magnitude of Eq 
in Eq. {EED. 



X. CONCLUSION 

A general implicit formula for the longitudinal kick due 
to the coherent synchrotron radiation has been developed 
for particles on a common orbit. This formalism will han- 
dle any geometry of bends and drifts. For simulations, 
this formula is to be preferred over the explicit ultra- 
relativistic formula because the implicit formula does not 
have a singularity at z = and is applicable at lower par- 
ticle energies and smaller length scales. 

Additionally, a heuristic formula for the longitudinal 
space charge kick has been presented which takes into 
account transverse displacements of the kicked particles. 
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FIG. 17: Numerical results using the A&Y code with varying 
dimensions of a rectangular chamber. A lm, 7.5 degree bend 
with a 2ps long bunch of 0.82nC charge was used. Contours 
represent the energy change induced due to the CSR wake 
field in 2000eV increments. The shielding effect thus not only 
due to the vertical dimension, but primarily due to the smaller 
of the two beam-pipe dimensions. 



It initially might seem surprising that vacuum cham- 
bers of many centimeter width can have a shielding effect 
on synchrotron radiation with much smaller wavelength. 
We therefore add a heuristic explanation here. 
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Starting with Schwinger [l|, shielding by a vacuum 
chamber has often be considered by studying infinite hor- 
izontal, conducting plates. But the following heuristic 
argument indicates why it is both the vertical and the 
horizontal boundary that determines shielding of CSR. 
Numerically solving Maxwell's equation in the vacuum 
chamber supports this argument as shown in Fig. |T7j 

A highly relativistic particle emits synchrotron radia- 
tion within a narrow cone. For the radiation's compo- 
nent of wavelength A, the opening angle of this cone is 
approximately AO = (7^)3 in the horizontal and vertical. 
The opening angle in the vertical determines the vertical 
width of the radiation load on the vacuum chamber wall, 
and the horizontal divergence can be observed by shining 
through a pinhole. 

The radiation field builds up within a radiation buildup 
time of At = ^ (^) ! - During this time, the radiation 
fields produced by the electron coherently add up to form 
the full radiation power. If the radiation does not in- 
terfere with an obstacle, for example the vacuum pipe, 



within this time, the electron looses as much energy as it 
would without any vacuum pipe. 

The width and hight of the radiation cone that builds 
up during the radiation buildup time is therefore given 
by w r ~ h r ~ cAtAO = R(^)^ . Vacuum chambers 
that have smaller dimensions interfere with the radiation 
process and shield the part of radiation for which 



A > Min 




(Al) 



Wavelengths are therefore shielded when they are above 
a length that is much smaller than the chamber dimen- 
sions. 

While we have used a very approximate heuristic argu- 
ment, Fig. [T71 computed by the A&Y code, indeed shows 
that both dimensions can lead to shielding, and that to 
first approximation only the smaller of the two dimen- 
sions is relevant. 
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